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Abstract 

We show that the degenerate sector of Spin (4) Plebanski formulation of four-dimensional 
gravity is exactly solvable and describes covariantly embedded SU(2) BF theory. This fact 
provides its spin foam quantization and allows to test various approaches of imposing the 
simplicity constraints. Our analysis suggests a unique method of imposing the constraints 
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which leads to a consistent and well defined spin foam model. 
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1 Introduction 

The path integral quantization is a powerful method of quantizing dynamical systems. However, 
in the case of gravity the usual path integral is not well-defined due to the well known non- 
renormalizability problem. In this situation the only hope is to use some non-perturbative 
methods. A particularly interesting possibility to construct a well-defined path integral for 
quantum gravity is suggested by the spin foam approach [HE]. In its framework, transition 
amplitudes between quantum states of spatial geometry are represented as sums over spin foams, 
i.e. two-dimensional cellular complexes colored by ceratin group theoretic data. 

In three dimensions this approach is realized by two spin foam models: Ponzano-Regge 
model [3] and Turaev-Viro model [I], which provide quantization for three-dimensional gravity 
with a vanishing and a positive cosmological constant, respectively. These models are well- 
grounded since they can be connected either to the loop quantization [5J or to the Chern-Simons 
quantization [HI |7], which has a solid basis based on the canonical analysis. Thus, in the three- 
dimensional case the spin foam quantization is the full-fledged approach to quantum gravity. 

In four dimensions the status of the spin foam quantization is much more controversial. 
The only models which are affordable to the direct spin foam quantization are the so-called BF 



1 



theories [5] which are described by the following action 

S BF [u,B] = [ d 4 xTr (B A F(u)) , (1.1) 

J M 

where u is the connection one- form valued in the Lie algebra q of a certain group G, F(u) is 
its curvature two-form, B is a two-form also valued in q, and Tr is evaluated using the Killing 
form on g. These theories are closed cousins of three-dimensional gravity. First of all, the action 
of the latter has a form similar to (11. ip with the only difference that B should be taken as a 
one-form and interpreted as a triad. And second, all these theories are actually topological. As a 
result, as in the case of three-dimensional gravity, their spin foam quantization can be obtained 
on a fixed triangulation on which it does not depend [9] (for recent developments, see also [TO]). 
In contrast, for non-topological theories, the dependence on the triangulation should be removed 
by summing over all of them or over their dual two-complexes. 

Besides these topological spin foam models of BF theories, there have been suggested several 
models supposed to describe gravity. Among them the most prominent ones are the Barrett- 
Crane (BC) model [HI [12] and the so-called new models, which are due to Engle, Pereira, Rovelli, 
Livine (EPRL) [13] and Freidel, Krasnov (FK) [TH [15]. However, all of them have been derived 
using certain assumptions and simplifications. Their common starting point is the so-called 
Plebanski theory, which represents gravity as a constrained BF theory with G = Spin (4) in the 
Riemannian or G = SL(2, C) in the Lorentzian case. Namely, it realizes the simple fact that, if 
the action (ll.ip is supplemented by a constraint ensuring that 

B=*(eAe), (1.2) 

it reproduces the usual Hilbert-Palatini formulation. The idea leading to the spin foam models 
mentioned above is that these constraints, called commonly simplicity, can be incorporated at 
the quantum level. Thus, they are supposed to be imposed on the spin foam representation of 
the quantum BF theory and should convert the trivial dynamics of a topological theory into that 
of quantum gravity. This strategy, which can be summarized as "first quantize, then constrain", 
is now the usual approach to four- dimensional spin foam models, and what distinguishes various 
models is only the way the simplicity constraints are incorporated. 

However, although this strategy seems to be well motivated and leads to interesting results, 
it does not agree with the rules of quantization of constraint systems. The simplicity constraints 
are known to be second class and therefore affect the symplectic structure to be quantized, which 
can be evaluated using the Dirac bracket. It has been argued that this and some other effects are 
not taken properly into account in the spin foam approach based on the above strategy [T6l [TT] . 
In particular, relying on the consistency with the canonical quantization, it has been suggested 
a certain modification of the vertex amplitude [18j, which is the most important quantity in spin 
foam models encoding their dynamics. If one follows the usual strategy, the vertex coincides 
with the one of BF theory, but restricted to a set of representations and intertwiners, assigned 
to the elements of the cellular decomposition, satisfying the simplicity constraints. Equivalently, 
it can be represented as the boundary state associated with a four-simplex evaluated on a flat 
connection, or as 
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where ^[gf\ is the boundary state depending on the group elements assigned to triangles of the 
four-simplex, the product goes over its tetrahedra, u(f)/d(f) denotes the upper/down tetra- 
hedron sharing triangle /, and the integration measure should be taken to be the usual Haar 
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measure on the group, T>g = dg. In [IS] it has been argued that the correct vertex for the 
constrained theory should be given by the same formula, but with the measure which involves 
the delta-function of secondary second class constraints, conjugate to the simplicity and ensuring 
that they are of second class. 

Given this situation when there are several proposals for the spin foam quantization of 
quantum gravity, it is desirable to have some simplified models which allow to test various 
features of these proposals. Moreover, since most problems and ambiguities arising in four- 
dimensional spin foam models come from the difficulties in imposing the simplicity constraints, 
such a model should mimic the structure of Plebanski formulation. In other words, it should be 
of the following form 

5 < Th2 = ^Thi + constraints, (1.4) 

where the constraints convert the theory Thl into the theory Th2. Finally, to be a good test- 
ground, both theories, given by the actions Srhi and Srh2, should have known spin foam rep- 
resentations. Then we can verify which of the methods to impose the constraints reduces the 
spin foam quantization of Thl to that of Th2. If a method does not allow to recover the known 
quantization of Th2, this strongly suggests that it is not applicable also in the case of Plebanski 
theory. On the other hand, if we find a quantization procedure which passes through our test, 
one can hope that it will work for gravity as well. 

In four dimensions the above requirements suggest that Thl and Th2 should be of BF type 
(11. ip because these are the only theories with a well-established spin foam quantization. Then 
the constraints can be used to reduce a gauge group G to its subgroup H C G. Thus, one arrives 
at an essentially unique model suitable for all our purposes: it should represent BF theory with 
the gauge group, for example, Spin(4) reduced by means of some constraints to the SU(2) BF 
theory. 

In [18] it has been noticed that such a reduction does take place if the -B-field and the spin 
connection are required to satisfy 

B IJ xj = } J\I l {x)uj kl = 2x [I dx J \ (1.5) 

where x 1 is a normal vector describing the embedding of SU(2) into Spin (4) and Jj/l( x ) * s 
the projector on the orthogonal completion of the su(2) subalgebra stabilizing x 1 (see appendix 
lA.lj) . Moreover, it was demonstrated that one recovers the known spin foam representation of the 
SU(2) BF theory from the known spin foam quantization of the Spin (4) BF theory only if one 
incorporates the constraints on the spin connection into the definition of the vertex amplitude 
in the way described below (II .3p . However, the important drawback of this consideration was 
that the constraints (ll.5p have been imposed by hand and not derived from a classical action of 
Plebanski type. 

An analogous model, but based on a solid canonical analysis, was proposed recently in [19] 
(see also [20]). It also represents the reduction of Spin(4) BF to SU(2) BF, but this time 
in three dimensions where the resulting theory is nothing else but three-dimensional gravity 
with vanishing cosmological constant. The constraints performing this reduction are a direct 
generalization of the simplicity constraints of the Plebanski formulation so that one could test 
whether one of the standard methods of imposing the simplicity reproduces the Ponzano-Regge 
model. It was shown that, as in four dimensions |18] . it becomes possible to recover the Ponzano- 
Regge vertex amplitude only if the simplicity constraints are supplemented by the secondary 
second class constraints restricting the holonomies of the spin connection. Furthermore, in [20J , 
the same model was reanalyzed and it was shown how it can be quantized so that the result 
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is consistent with the canonical quantization approach and reproduces all ingredients of the 
Ponzano-Regge model, including face, edge and vertex amplitudes. 

Although these results are already very suggestive, they may be argued to be not directly 
applicable to the four- dimensional spin foam models. The reason for this is that in three di- 
mensions, for spin foams dual to a simplicial discretization, all edges are three- valent. As a 
result, the space of gauge invariant intertwiners assigned to these edges is one- dimensional and 
the most interesting part of the simplicity constraints, the so-called cross simplicity, is automat- 
ically satisfied. On the other hand, it is the treatment of the cross simplicity that distinguishes 
the spin foam models mentioned above, so that one could think that this part of the constraints 
is particularly crucial for the whole construction. 

To avoid such critics, in this paper we return to the original four-dimensional model of [IS] 
and fill there all the gaps both at classical and quantum level. Namely, in section [2] we propose 
a classical action which represents the Spin(4) BF theory reduced down to SU(2) BF, where the 
constraints (jl.5p appear as primary and secondary constraints, respectively. Remarkably, this 
action represents a system very close to the physical system we are interested in — it describes 
the degenerate sector of Plebanski theory. Thus, the latter is exactly solvable and can be seen as 
a covariant embedding of the well known topological theory. We provide a thorough canonical 
analysis of this degenerate sector both with a partial gauge fixing and without it. In particular, 
we show that adding the degeneracy condition to the usual simplicity constraints, one generates 
constraints for constraints. They are responsible for the well known fact [2T1 [22] that the phase 
space of degenerate configurations is larger than its non-degenerate version. It is also possible 
to include the Immirzi parameter into the model which does not lead to any complications. 

Then in section [3] we consider the spin foam quantization of this model. First, we apply 
the usual quantization strategy employed in the EPRL and FK approaches. Since our classical 
action differs from Plebanski theory only by the presence of the degeneracy condition, to get 
its spin foam quantization, it is sufficient to extract the degenerate sector of the new spin foam 
models. If it is done, following the usual ideas, by restricting the boundary or kinematical data, 
the result strongly disagrees with the known vertex amplitude for the four- dimensional SU(2) 
BF theory. On the other hand, if the degeneracy condition is represented as a constraint on 
the group elements appearing in the integral formula (11. 3p for the vertex and is added to the 
integration measure, one does get the right result. Since this constraint can be equally seen as 
a discretization of the secondary second class constraints, this modification of the measure is in 
the full agreement with our proposal for the vertex amplitude spelled above. 

At the same time, the analysis of the constraint imposition in the framework of the new 
models shows that the constraint on holonomies ensuring the correct vertex amplitude erases 
all information about the solution of the original simplicity constraints of the EPRL and FK 
models. Thus, the main ingredients of these models appear to be irrelevant for getting the right 
dynamics, which calls for a reconsideration of these approaches. As an alternative, we suggest 
another quantization procedure, which summarizes the analysis of [18j [20] and is consistent with 
the canonical approach by construction. Being applied to the model under consideration, it gives 
precisely the right result: the Crane- Yetter model with the SU{2) structure group. Moreover, 
in the course of evaluation of the partition function, we clarify the role of different constraints 
in the spin foam quantization. In particular, we observe that the vertex amplitude is completely 
determined by the secondary second class constraints putting restrictions on holonomies of the 
spin connection, whereas the primary simplicity constraints affect only the gluing of different 
vertex contributions and are not relevant for dynamics. A discussion of these and other issues 
can be found in the concluding section. 
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Our conventions are explained in appendix |A] We restrict ourselves to the Riemannian case 
to not bother the reader with signs which otherwise would pop out here and there. However, 
all the presented results are easily generalized to the Lorentzian case as well. In appendix IB1 we 
present the details of the canonical analysis of the degenerate Plebanski sector without a partial 
gauge fixing, whereas the last appendix O makes explicit the constraints for constraints. 

2 Degenerate sector of Plebanski formulation 

2.1 The action and constraints 

Let us start from the usual Plebanski action 

S dPl [Lu,B,\] = l [ d 4 x 

where the Lagrange multiplier field A is chosen to be a spacetime pseudo-tensor satisfying 
the following symmetry properties: it is antisymmetric in the first and second pair of indices, 
y^vpa _ ^[fj,u][ P a]^ an( j - g S y mme tric under their exchange, [fxu] -H- [pa]. Usually, one also adds the 
tracelessness condition e^ vpcT \^ upa = which we however omit. As a result, the variation with 
respect to the Lagrange multiplier generates the following constraints 

$iiup<r = ^ £ UKLB I J / Bf a L = 0. (2.2) 

They represent the usual 20 simplicity constraints supplemented by an additional condition which 
forces the S-field to belong to the degenerate sector, i.e. to give a vanishing four-dimensional 
volume V = j% e pupa ^^ upa . However, in contrast to the usual case, not all of these 21 constraints 
are independent. It turns out that there are 6 constraints for constraints, which we display 
explicitly in appendix O As a result, one remains only with 15 independent constraints which 
thus reduce the number of independent components of the i?-field from 36 to 21. 

It is easy to see that there are two sectors of solutions to the simplicity constrains ( 12. 2D . 

"deg-gravitational" : Bj£ = -e IJ KL x K b L ^ (2.3a) 
"deg-topological" : B% = xH$,, (2.3b) 

where the vector x 1 is supposed to be normalized as x 1 xj = 1. In analogy with the usual case, 
we call them "degenerate gravitational" and "degenerate topological" Q The reason for this will 
be clear from what follows. Both solutions (12. 3D contain 21 independent degrees of freedom, i.e. 
the same as the solution space of the simplicity constraints: 3 degrees of freedom are contained 
in x 1 and 18 are described by b 1 because the latter field can be chosen to satisfy the linear 
constraints 

xjb 1 ^ = 0. (2.4) 

They fix uniquely the ambiguity in b 1 and will be always assumed to hold in the following 
analysis. 

1 There might be also a "twice degenerate" sector because the proof in appendix [C] that there are only 6 
constraints for constraints relies on the assumption that a certain matrix constructed from the -B-field is invertible. 
But we do not consider here this possibility and assume that the field b* is generic so that it ensures certain 
non-degeneracy conditions appearing in the course of our analysis. 
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If one fixes the vector x 1 , thereby reducing the gauge symmetry from Spin(4) to the sub- 
group SU(2) X which preserves this vector, the two sectors of solutions (12.31) can be equivalently 
characterized by linear simplicity constraints 

"deg-gravitational" : $j^ )J = xjBf£ = 0, (2.5a) 
"deg-topological" : = e IJ KL xjB% L = 0, (2.5b) 

again in direct analogy with the usual non-degenerate case [23], [2U [TU [25]|^| However, the 
difference is that here the linear constraints (12.51) exhaust all simplicity constraints, whereas in 
the non-degenerate case they should be supplemented by the volume constraint. Indeed, since 
both constraints (12. 5p satisfy xi&^y 1 = x/$^° p ^ = 0, they give rise to 18 independent equations. 
This reduces the number of independent components of the -B-field to 18, which coincides with 
the number of independent degrees of freedom described by b^. 

To understand the meaning of the two sectors, let us assume for simplicity that x 1 = const. 
Then one can easily extract an equation without derivatives from the equation of motion obtained 
by variation of (12 .ip with respect to the spin connection 

£ w DvB u = o, (2.6) 

where _D M is the covariant derivative defined by u^. Indeed, let us contract the equation (12. 6ft 
with xj in the deg-gravitational sector and with e KL ijXl in the deg-topological sector. Taking 
into account the form of the B- field (12. 3p . in both cases one finds 

e^e IJ KL xjb%{^ M x M ) = 0. (2.7) 

Under the assumption of invertibility of the matrix multiplying the spin connection, or more 
precisely xjuj/ , this equation requires that 

xjul J = =► xjF^ = 0. (2.8) 

This result indicates that only the part of the connection describing the SU(2)a; subgroup sur- 
vives, whereas the orthogonal part vanishes. Now we can plug the solution of the simplicity 
constraints (12.31) into the original action (12.11) . It is immediate to see that, due to (12. 8p . in the 
deg-topological sector the resulting action identically vanishes, whereas in the deg-gravitational 
sector it becomes 

S v [u,b] = \f d 4 i^%^i^. (2-9) 
4 Jm 

This is nothing else but the action of the four-dimensional SU(2) BF theory covariantly embedded 
into Spin(4). Fixing the time gauge x 1 = 5$, one recovers its usual form written in terms of 
the 5u(2)-valued 2-form b % ^ v and the su(2)-connection u^. The other components of the original 
fields vanish due to O and (jUg) . 

Since in the deg-topological sector we do not obtain any meaningful theory, we will be mainly 
concentrated on the deg-gravitational sector. Remarkably, it is given by a well known topological 
theory so that we know its exact classical as well as quantum description. For our purposes it 
will be however important to understand also the canonical structure of the original Spin(4)- 
invariant theory ( 12. ip which we present in the next subsection. Then in subsection 12.31 we will 
see how the reduction to the SU(2) BF theory described here is established in the Hamiltonian 
formalism. The reader who is not interested in details of this canonical analysis can proceed 
directly to section 12.41 

2 Note however that in the non-degenerate case the roles of $( gr ) and $( to P) are exchanged, i.e. the former 
corresponds to the topological sector and the latter to the gravitational. 
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2.2 Canonical analysis 



Here we present the Hamiltonian formulation of the degenerate gravitational sector of Plebanski 
theory in a partially fixed gauge. Namely, we fix the normal field x I (x) to be a given function 
on spacetime. The choice of the normal allows to replace the quadratic simplicity constraints 
(12.21) by their linearized version (I2.5p . which can be done directly in the action and gives the 
possibility to restrict to the particular sector we are interested in from the very beginning. An 
analogous formulation of the gravitational sector of Plebanski theory has been considered in 
[25] (see also [26]). However, the important difference of our approach is that the normal x 1 is 
considered as a fixed non- dynamical variable. This is motivated by the following application of 
these results to the spin foam quantization of our model. This quantization is implemented via 
a path integral where the gauge freedom generated by boost transformations should be fixed by 
a gauge choice. The most convenient way to do this is precisely to fix the normal x 1 . This gauge 
fixing is analogous to what is done in the standard loop quantum gravity where one imposes the 
so-called time gauge corresponding to a particular choice of x 1 = 5$. Here we could also restrict 
ourselves to this simple gauge, in which case the following derivation considerably simplifies. We 
however prefer to keep x 1 an arbitrary function to show that in this general case one obtains 
nice covariant structures. 

We emphasize that the results presented here can also be derived going through the complete 
canonical analysis of the original action ( 12. ip carried out without imposing any gauge fixing. We 
provide such analysis in appendix [B] where the canonical structure of both solution sectors ( 12. 3 j) 
is elucidated. They possess a very intricate constraint structure which however reduces to the 
one of this subsection upon restricting to the deg-gravitational sector, fixing x 1 , and solving 
auxiliary constraints. 

Thus, our starting point is the following action 

S deg - gr k B,\-x] = \f d 4 x [e^B^FH + iXfxjBffl , (2.10) 

where x 1 is just a parameter and not a dynamical variable. After the 3+1 decomposition of this 
action, one can recognize that the phase space is parametrized by 

ui J and P? J = e abc B I b J (2.11) 

with the canonical commutation relations 

R /J (x),pL(y)} = <*aO(*,y)- (2-i2) 

The variables Uq J , B^ and appear without time derivatives and therefore play the role of 
Lagrange multipliers. However, not all of them generate constraints. Whereas the variation with 
respect to ul J and does give rise to the primary constraints 

Qu = D a P?j^0, (2.13a) 
= x J PjjfaQ, (2.13b) 

the variation with respect to A^ a leads to a condition on the Lagrange multipliers, 

XjBH = 0. (2.14) 

At the same time, B 1 ^ gives the following equation 

e ahc F^ + A\\?x A « 0. (2.15) 
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It can be split into two parts: one gives a condition on the Lagrange multipliers 

\T = -\s abc F^xj, (2.16) 
and the second is a primary constraint 

C/° = e IJKL x J e abc F b K c L « 0. (2.17) 

Taking into account the condition (12.141) . the Hamiltonian is given by a linear combination of 
the primary constraints introduced above 

-H= l - e IJ KL xjB« L C} + u^Qu + e abc Xf^. (2.18) 

Now one should find the conditions under which the primary constraints are preserved in 
time. Since their evolution is generated by the Hamiltonian (I2.18p . this boils down to the study 
of the constraint algebra. Introducing the smeared constraints 



J dWC?, 



Q(n) = J d^n^Qu, C a (v) = / dWC 7 a , (2.19) 

the non-trivial commutators are given by 

{S(n),G(m)} = Q([n,m]), 

{g(n),C a (v)} =C a (n-v)-A j ^e IJ KL \fxjv K n LM x M , 

{G(n), $?} = -n IJ ^j + Phn JK x Kl 
{C»,<^} =-2e IJKL v J x K e abc D c x L , 



(2.20) 



where we used (I2.16p . To apply these results, we start from the primary constraint Its 
conservation leads to the condition 

6? = 2e ahc B^D c xj + P?j {u J Q K x K + x J ) « 0. (2.21) 

Combined with the three equations following from the Gauss constraint 

x J Gu = D a § a j - P?jD a x J « 0, (2.22) 

it gives rise to the condition on the Lagrange multipliers uJq J 

Dqx 1 = (2.23) 

and to the secondary constraints 

^ a = DaX 1 « 0. (2.24) 

Then the conservation of the Gauss constraint Qu gives a relation between the Lagrange mul- 
tipliers X^ u 

g IJ ^-4x^B J ] a K X K a -e abc x [I P^ = =► 5^ = 0, (2.25) 
whereas the conservation of Cj does not generate new conditions. 
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The next step is to study the secondary constraints (12.241) . In fact, due to the identity 
xi^a = 0, they give only 9 independent equations. These constraints satisfy the following 
commutation relations 

¥ b } = 0, {G(n), ¥ b } = xjD b n IJ , * 6 J } = -~ - xjx J )5(x, y). (2.26) 

Taking into account the restriction (12.231) . the conservation of *fr a therefore amounts to vanishing 
of all Lagrange multipliers Xf = 0. (Recall that these multipliers can be chosen to satisfy 
x J A" b = from the very beginning so that the number of their independent components equals 
the number of independent secondary constraints.) Moreover, it is easy to check that 

£ *bc F u Xj = 2e abc D b V I c ==> A? a = 0, (2.27) 

which fixes the last undetermined part of the Lagrange multipliers A/". Thus, the stabilization 
procedure stops at this point. 

Once we found all the constraints, we can study whether they are of first or second class. 
Due to (I2.22p . we remain only with four types of constraints: 

Cf, g I = e IJ KL x J g KL , ¥ a . (2.28) 

Furthermore, due to the Bianchi identity, the first ones additionally satisfy 

D a Cf = WtCfaj] + 2e IJKL x J e abc ^D b m L c . (2.29) 

Using the commutation relations (I2.20p and ( I2.26p . it is trivial to check that the constraints 
Cf and Qi are first class, whereas and ^ are second class. As a result, the 18 + 18 = 36 
configuration variables are restricted by (9 — 3) + 3 = 9 first class and 9 + 9 = 18 second class 
constraints, which leaves us with a zero dimensional physical phase space. This confirms that 
this theory is topological, i.e. it does not contain propagating degrees of freedom. 

2.3 Reduction to SU(2) BF theory 

The meaning of the constraints (12.281) is quite transparent: the second class constraints, 
and \E^, fix the off-diagonal (boost) degrees of freedom in the chiral decomposition of the so (4) 
algebra or, more precisely, the configuration variables from the orthogonal completion to the 
su(2) x subalgebra. The remaining constraints, which describe the dynamics of the variables 
from this subalgebra, are nothing else but the usual constraints of the SU(2) BF theory |27J. 
Namely, Cf gives the flatness condition for an SU(2)-connection and Qi is the corresponding 
Gauss constraint generating SU(2) gauge transformations. This becomes especially clear in the 
time gauge x 1 = 5$ where the second class constraints simply mean that P^ a = P^ a and 
ioi + ^ = OJa . For a constant x 1 these relations get rotated by a constant Spin(4) transformation 
mapping 5q into x 1 . It is however amusing to see how they generalize to the case of arbitrary 
gauge where x 1 can vary in spacetime. 

Let us introduce the projections of our fields on the su(2) x subalgebra 

Ph = IuWPkl, b£ = I&WBS?, (2.30) 

where lfj L (x) is the projector given in (1A.2I) . so that the new fields solve (12.13bl) and (12.141) . 
respectively. Besides, we define the following connection 

K J = I kl{xW» +2x^d,x J \ (2.31) 
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where the last term takes care about variations of the normal field. This connection coincides 
with the original spin connection uf/ on the surface of (12.231) and (I2.24p and satisfies the con- 
straint 

J I K J L^ L = 2^ I d,x J \ (2.32) 

which is identical to the one appearing in the Lorentz covariant formulation of loop quantum 
gravity [2"8ll2"U] . The characteristic feature of this connection is that its holonomies map a vector 
from su(2) Xl to a vector in su(2) X2 subalgebra [16] , i.e. for constant x 1 they belong to the SU(2) a . 
subgroup. 

In terms of the new variables and taking into account all conditions on the Lagrange multi- 
pliers, the 3+1 decomposed action ( 12.101) can be written as 



Sdeg-gr = I tit I d 3 x 



(2.33) 



where F^(A) is the curvature of the connection A 1 ^ and 

Gi? = d a p a IJ + [A a ,p a ) IJ . (2.34) 

Setting the second class constraints to zero, one obtains the usual BF action where all variables 
are projected down to SU(2) X . It provides a covariant embedding of the SU(2) BF theory into 
the Spin (4) formalism. 



2.4 Summary 

Let us summarize what we have found. The theory (12.11) describing the degenerate sector of Ple- 
banski formulation of general relativity has two sectors of solutions of the simplicity constraints 
( 12. 2ft . In the sector which we called "degenerate gravitational", it reduces to the four-dimensional 
SU(2) BF theory covariantly embedded into Spin(4) gauge group. The embedding is character- 
ized by the normal vector x 1 . In the partially fixed gauge where x 1 is a fixed function, the theory 
possesses two types of second class constraints conjugate to each other: 

primary = x j Pjj, secondary = D a x J . (2.35) 

The remaining constraints are first class and generate the gauge symmetries of the SU(2) BF 
theory. 

The presence of the second class constraints, as usual, leads to a modification of the symplectic 
structure: the Poisson bracket (12. 12ft has to be replaced by the appropriate Dirac bracket. The 
latter can be easily calculated and is given by 

U J (x),P b KL(y)}n = <M J L (x)<Kx,y) (2-36) 

with all other commutators being vanishing. This result again demonstrates that only the SU(2) X 
part of the configuration variables is dynamical and implies that Pfj and uJ a J are not canonically 
conjugate anymore. 
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2.5 Inclusion of the Immirzi parameter 



It is easy to include the Immirzi parameter [30] into our model. To this end, one makes the 
usual replacement of the I?-field in the BF part of the action by the combination B + - * B. It 
leads to a mixing of the two solution sectors (12. 3p . Now both of them reduce to the SU(2) BF 
theory, just in the deg-topological case the resulting action is multiplied by the factor I/7. 
In the Hamiltonian formulation of the deg-gravitational sector described by the action 



0(7) 
o. 



deg-gr 



[u, B, A; x] 



d 4 x 



M 



tdIJ , 1 A J n KL 



Fit + A\fxjB IJ 



per 



(IV 



(2.37) 



the presence of the Immirzi parameter affects the Poisson symplectic structure so that the new 
canonical variables are 



to 



1 j 



and 



pa 



(1 + Pfj. 



(2.38) 



Nevertheless, all the remaining structure does not change and there are just slight modifications 
in the stabilization procedure of section 12.21 In particular, the final constraints acting on the 
phase space (12.381) comprise Cf and Qj which are first class and which are second class, 

and all these constraints are given by the same expressions as without 7. Moreover, due to 
the second class constraints, the 7-dependent symplectic structure given by Poisson brackets is 
replaced by the 7-independent one described by the same Dirac brackets (12.361) as before. 



3 Spin foam quantization 

In this section we consider the spin foam quantization of our model. Since after implementing 
the second class constraints (in the deg-gravitational sector) it coincides with the SU(2) BF 
theory, we know what the final result should be: it is given by the Crane- Yetter model [3"Tl E] 
with the structure group SU(2) represented by the following spin foam state sum 

^ Y (2) (A*) = e e n &f + 1) n ^ su(2) > 

i^e feA* v£A* 

where A* is a two-complex dual to a simplicial triangulation A of the spacetime manifold A4, 
jf labels irreducible representations of SU(2) attached to the faces of A*, i e are SU(2) invariant 
intertwiners assigned to its edges, and A v is the vertex amplitude given by the SU(2) {15j} 
symbol. The latter is obtained by evaluation of the SU(2) spin network represented by the 
pentagon graph, which is dual to the boundary of a 4-simplex a G A dual to the vertex v £ A*, 

h 

A s v ui2) m = im = , / J24 \. — > ^ (3-2) 

\\ji:t/ / 

\ /^ 25 ^ x y\ / 
Ji2\ / \<r \ A' 4 5 



11 J51 ^ 
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The normalization of the intertwiners used in this evaluation is defined in appendix IA.21 

However, we would like to proceed in a different way which would avoid quantizing the degrees 
of freedom on the constraint surface only. Our aim is to find a quantization of the original theory 
with constraints, (12. ip or f 1 2 . 1 j) (or even (I2.37P ). such that it reproduces the SU(2) Crane- Yetter 
model (13. ip . In particular, we would like to check whether the quantization strategies used to 
get the EPRL or the FK model are able to do this. Since our model is only slightly different 
from the gravitational sector of Plebanski theory, our study represents a very serious test on the 
validity of these quantization approaches. 

3.1 Discretization 

Before we start discussing the quantization, let us put our model on a cellular complex and 
discretize its variables. As usual, the -B-field is discretized by associating Lie algebra valued 
elements Bf e so (4) to the faces of the dual complex A* and can be obtained as integrals of the 
5-field over the dual triangles 



The spin connection u gives rise to group elements g e which coincide with its holonomies along 
the edges e G A*. However, it is also convenient to introduce the holonomies g ve , going from 
vertex v to the center of edge e, and gj e , going from the center of face / also to the center of 
the edge. g ev and g e j will denote their inverse. The former provide a refined version of our basic 
dynamical variables which are obtained as g e = g ve g e v', where v and v' are the two vertices joined 
by e. On the other hand, gf e are needed to bring the bivectors (13. 3p to the reference frame where 
g ve are acting. In particular, we define 



In this sense, these group elements may be considered as non-dynamical auxiliary variables com- 
pleting the definition of the discrete 5-field. Altogether, g ve and gf e provide the discretization 
of the spin connection on the two-complex obtained by subdivision of A* into wedges, which are 
in one-to-one correspondence with pairs (vf), as shown on Fig. [1] [32J. 

Finally, we should discretize the normal field x 1 appearing explicitly in the gauge fixed 
action (12.101) and in the linearized simplicity constraints (12. 5ft . This field is analogous to a 
similar field appearing in the canonical formulation of the usual Hoist and Plebanski actions [33] 
where it describes the normal to the three-dimensional spacelike slices. In spin foam models it 
usually gives rise to the normal vectors x e , which can be viewed as elements of the factor space 
X = Spin(4)/SU(2) [231 HEl ESI EI] • Such a vector is interpreted geometrically as the normal to 
the tetrahedron dual to edge e. However, in the degenerate case it is more natural to associate 
such normal vectors to 4-simplices. Thus, at the discrete level the field x 1 will be represented 
by a set of unit vectors x v . 

With these definitions we can now discretize the simplicity constraints. Their quadratic 
form f]2. 21) is discretized as usual giving rise to diagonal, cross and volume simplicity obtained 
by averaging the two bivectors over the same triangle, different triangles belonging to the same 
tetrahedron, or non-intersecting triangles of the same 4-simplex, respectively. The only difference 
with the usual case is the form of the volume constraint which requires that the geometric volume 
of the 4-simplex vanishes. 

On the other hand, the discretization of the linear simplicity constraints (I2.5al) involves an 
extra ingredient. Indeed, they relate the fields, x 1 and B IJ , which after discretization live at 




(3-3) 



B e f = g e fB f g fe . 



(3-4) 
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Figure 1: Two faces sharing edge e! and divided into wedges. The small arrow indicates the 
orientation of the face and the big arrow shows how one evaluates the curvature associated to 
the wedge. For a given vertex v, the two adjacent edges e and e' are denoted as u(f) and d(f), 
according to the orientation of / to which they both belong. 

different elements of the cellular decomposition: at vertices and faces, respectively. Due to this, 
the bivectors should be transported to the reference frame of a vertex using the holonomies 
introduced above. As a result, the discrete simplicity constraints read as follows 

g ve B ef g ev ■ x v = 0, for V/DeDu. (3.5) 

It is important to notice that since there are two ways to connect a face to a vertex (by going 
either through e = u(f) or e = d(f)), the conditions (I3.5p constrain not only the bivectors, but 
also the holonomies. This shows that at the discrete level the primary and secondary constraints 
are not well distinguished from each other. 

To make contact with the new spin foam models, it is convenient to change a bit the point 
of view and to write the simplicity constraints in the reference frame of a tetrahedron or its dual 
edge. To this end, we define 

Xe{v) = QevX-u (3.6) 

so that the condition (13.51) becomes 

Blj{x e {v))j = 0, for V/DeDt). (3.7) 

Up to insertion of the Hodge operator, this is the usual form of the linear simplicity constraints 
used in the new spin foam models. In our case it should be supplemented by the additional 
requirement that the normals x e (v) originate from the same vector x v and therefore must satisfy 

g ve x e (v) = g ve >x e >(v), for Ve, e' D v. (3.8) 

Then we turn to the secondary second class constraints (12.24)) . They restrict the holonomies 
of the spin connection and at the discrete level read as 

x v = g e x' v , for e D v, v' (3.9) 
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Being combined with (I3.6p . they can be equivalently rewritten as 



x e (v) = x e (v'). (3.10) 

This shows that, provided the secondary constraints are imposed, one can drop the v -dependence 
of the normals x e . Moreover, this suggests to introduce the normal vectors associated to all 
elements of the spin foam cellular complex: x v , x e and Xf. This pluralistic point of view allows 
to formulate all constraints in a simple and uniform way. Indeed, they become equivalent to the 
following relations 

B I f J (x f ) J = 0, (3.11a) 
Xe = 9evXv, x f = g fe x e . (3.11b) 

These relations represent the straightforward discretization of the primary and secondary second 
class constraints (I2.35p . respectively. We emphasize that it is crucial to consider the primary 
and secondary constraints on equal footing. For example, taken alone, ( 13. Haft is not sufficient to 
generate the degeneracy condition for a 4-simplex. On the other hand, altogether the conditions 
( 13. lip provide an elegant and amazingly simple discrete formulation of all second class constraints 
of the continuous theory. 

To facilitate the use of the constraints in the discretized path integral, let us rewrite them 
using the chiral decomposition of Spin(4) (see appendix IA.1I) . To this end, we note that each 
normal vector i6l gives rise to an element x G SU(2) defined by 

x = g-(gir\ (3-12) 

where g x G Spin(4) is a representative of x which choice does not affect the definition of x. Then 
it is clear that ( 13. lip is equivalent to 

Bj=x f B+x]\ (3.13a) 
9ve = x^+xj 1 , gj e = x f gj e x-\ (3.13b) 

The other constraints appearing above like (13. 7p and (13.81) have a similar representation. 
3.2 The usual strategy 

Our first aim is to apply to our model the quantization strategy employed in the EPRL and FK 
spin foam models and based on the idea "first quantize, then constrain" . This implies that one 
should start from the unconstrained Spin (4) BF theory and incorporate the simplicity constraints 
at quantum level. The spin foam quantization of the unconstrained theory is provided by the 
Crane- Yetter model [9] with the structure group Spin(4) represented by the following state sum 

2c p Y n(4) (A*) = £ EII Ml + Wff + !) II {w+Hun, (3.14) 

\=(j+ ,j-)-+f Z=(i+ ,i-)-*e feA* veA* 

where the sum goes over all Spin(4) irreducible representations A/ and all Spin(4) invariant 
intertwiners X e . The simplicity constraints in this approach are supposed to restrict the allowed 
set of representations and intertwiners such that the resulting state sum provides the discretized 
path integral for the constrained theory. In our case, the result must reproduce the partition 
function (I3.ip . 
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However, it is clear that if the effect of the constraint imposition is only the reduction of the 
admissible group theoretic data, the state sum f l3.14|) will never reduce to (13. ip . Indeed, the two 
partition functions have different vertex amplitudes and it is impossible to reduce one to another 
by restricting the kinematical data. 

This approach can be realized explicitly by proceeding as follows. As is well known, the 
asymptotic analysis of both EPRL and FK models reveals that they contain a degenerate sector 
|22[ 135] . Its geometric interpretation precisely corresponds to the classical geometries described 
by our model. Thus, a simple way to get a spin foam quantization of ( I2.37P is to extract 
the degenerate sector from the EPRL or FK state sum. This can be achieved by expressing 
this state sum as an integral over coherent states [T4"l [T5] . It has been shown [35] that each 
quasiclassical Regge geometry contributing to the asymptotics of this integral can be uniquely 
reconstructed from the boundary data consisting of the SU(2) representations associated to 
triangles of A (or faces of A*) and coherent states assigned to each pair (ef). The degenerate 
sector is then extracted by restricting only to those boundary data which lead to degenerate 
Regge geometries. Since this restriction is imposed only on the boundary data, it affects only 
the set of representations and intertwiners which survive in the partition function (I3.14p after 
imposing the simplicity constraints. It does not affect the general form of the vertex amplitude 
and therefore cannot reproduce the desired result (13.21) . 

This method of imposing of the degeneracy condition has a clear drawback: it ensures the 
vanishing of the geometric volume of 4-simplices only in the quasiclassical limit. On the other 
hand, this condition appears as a part of our simplicity constraints which are expected to hold 
at the full quantum level. This suggests that one should look for an alternative approach. Such 
an approach does exist and moreover it leads to the correct vertex amplitude. But to realize it, 
we should be ready to go beyond the usual strategy and to constrain not only the kinematical 
data in (I3.14p . but also the group elements entering the definition of the vertex amplitude. 

In the previous section we showed that the simplicity constraints can be represented as 
a combination of two conditions, (13. 7p and (13. 8p . The former are the usual linear simplicity 
constraints in the non-degenerate sector of Plebanski formulation. The new spin foam models 
provide a way to implement them at quantum level so that this step can be considered as 
being already accomplished. Thus, it remains only to incorporate the second condition, which 
restricts us precisely to the degenerate sector by requiring that the normals to all tetrahedra of 
a 4-simplex, transferred to the same frame, coincide. 

Rather than a constraint on the normals x e (v), the condition (13. 8 p can be viewed as a 
constraint on the holonomies g ve . These holonomies coincide with the group elements appearing 
in the integral formula for the vertex amplitude of the new models 

A { jH] f ,k ef ,i e ) = / n d ^ e<S (r.,A(T)( jV ), fce/ ,i e ) [9vu(f)9 vd (f)' x e( v )] > ( 3 - 15 ) 

eDv 

where u(f) and d(f) denote the two edges belonging to the face / and sharing the vertex v (one 
of them is considered as "up" and the other as "down", see Fig. [TJ, and «S(r„,A/,fc e /,i e ) [g/^e] 
is the so-called projected spin network [3B] defined on the graph r„ dual to the boundary of a 
4-simplex (this is the same graph as the one appearing in (I3.2p ). The projected spin network is 
labeled by Spin(4) representations A/ attached to the links of the graph, SU(2) representations 
k e f assigned to the ends of the links, and SU(2) invariant intertwiners i e associated to the nodes. 
In both EPRL and FK models, the imposition of the simplicity constraints (13.71) leads to that 
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the Spin(4) representations Xf are defined in terms of the SU(2) representations jf a^H 

A (7) (J)= (^(l + 7)j,^|l-7|j)- (3.16) 

On the other hand, the representations k e f are treated differently: in the FK model they can 
be arbitrary, whereas in the EPRL model they are fixed to be k e f = jf. (For 7 = 0, the EPRL 
prescription gives X^°'(j) = k e f = and reproduces the BC model.) 

Given the integral representation (13.151) . the natural idea to incorporate the condition (13. 8p 
is to insert it into the measure. This amounts to adding the factor S ((^~ e )~ 1 x„5f+x~ 1 (f)) and 
it is straightforward to evaluate the resulting vertex amplitude. To this end, let us recall the 
explicit expression for the projected spin network 

S<P^M fe/^e] = <g)ie • (g^'^'D^ (< /)g /^ (/) ) *«'>', (3.17) 

where T)^ Xf ^(g) is the image of g G Spin(4) in representation Xf = (jf,jj) and C nni2, is the 
invariant map whose matrix elements are given by the Clebsch-Gordan coefficients. Then using 
the following property of the SU(2) matrix elements 

E <w£ D ™ (*p£3 (h) ciit = D ^ co > ( 3 - 18 ) 

m,m' ,n,n' 

where dj — 2j + 1, one can easily show that 

A^(j f ,k ef ,i e ) = /n^dfc^^-^ff+x: 1 ^)) 

= {15i}n<W,c/)/> ( 3 - 19 ) 

where the {15j } symbol is constructed out of ten representations kf — = kd(f)f and five 

representations characterizing the intertwiners i e , exactly as in (13 .21) . Thus, up to a normalization 
factor, we reproduced the correct vertex amplitude of the SU(2) Crane- Yetter model (13. ip ! 

This beautiful result seems to indicate that the quantization strategy realized by the new spin 
foam models passes our consistency check. However, its close inspection raises several questions. 
First of all, a striking feature of (13.191) is that it does not depend at all on the restrictions on the 
representations obtained by imposing the first part of the simplicity constraints (13. 7ft : it holds 
for any set of A/ and k e f. On the other hand, these restrictions are at the core of the new spin 
foam models and it is very puzzling that, after imposing the remaining part of the simplicity 
given by (13.81) . all information about them is completely erased. The only case where some 
information remains corresponds to 7 = in the EPRL model. But it is even worse. In this 
case kf are fixed to be zero and we do not get the correct vertex at all. Furthermore, although 
the final result (13.191) is perfectly fine for all values of the Immirzi parameter, the intermediate 
step (13.151) is not defined for irrational 7. 

3 One should take into account that the constraints Q3.7P differ from the linear simplicity constraints describing 
the gravitational sector in the new spin foam models by the absence of the Hodge operator, i.e. they are actually 
analogous to the constraints specifying the topological sector (see footnote [21) ■ This difference can be accounted 
by replacing the Immirzi parameter by its inverse. 



s, 



(r v ,x(t)(j f ),k ef ,i e ) 



9vu{J)9vd(f)i X e\ V ) 
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All these issues demonstrate that the constraints put forward by the EPRL and FK models 
are not really relevant for getting the correct spin foam dynamics, and even generate some 
strange unphysical effects like the quantization of the Immirzi parameter. On the other hand, 
the correct dynamics is obtained by imposing the constraints missing in the usual approach. In 
fact, there is a crucial difference between (13. 7p and ( 13. 8ft : although they both needed to discretize 
the primary simplicity constraints, the latter are better seen as a discretization of the secondary 
constraints (I2.24p . (Let us recall that at the discrete level there is no a clear distinction between 
the two types of constraints such that exists in the continuum theory.) Thus, the insertion 
of these constraints into the integration measure is exactly what has been suggested in the 
introduction (see below ( 11. 3p ). in our previous works [HJ[T6], an d in [19]. There it was claimed 
that the measure over holonomies should include the delta function of the secondary second class 
constraints. Here we see quite explicitly that this insertion is indeed necessary. 

As a result, we arrive at the following situation. It is indeed possible to extract from the 
EPRL and FK models the correct dynamics of the degenerate sector provided we incorporate the 
secondary second class constraints directly into the definition of the vertex amplitude. However, 
this modification of the vertex makes the imposition of the primary simplicity constraints com- 
pletely irrelevant. This questions the ability of these approaches to capture the right dynamics 
in the gravitational sector where the secondary constraints have been ignored so far. 

In fact, as we will show in the next subsections, there is a consistent way of quantizing the 
theory ( 12.3TI) which leads to the full correct result ( I3.ip . including not only the vertex, but 
also the edge and face amplitudes. It requires to take into account all constraints ( 13.111) . As a 
consequence, one can use any version of the primary simplicity constraints: either at vertices 
(13.51) . or at edges (13. 7p . or even at faces (13.11 al) — they all become equivalent. However, these 
constraints become important only for gluing the contributions of different simplices, whereas 
the vertex amplitude associated with a given simplex turns out to be completely determined 
by the secondary constraints restricting holonomies. These results clearly show that the usual 
strategy to the spin foam quantization, based on the use of only the primary constraints, is not 
satisfactory. 

3.3 Canonically inspired quantization 

In this subsection we provide the rules to construct the partition function of a constrained theory 
of Plebanski type. These rules summarize the results obtained for the vertex amplitude in [18] 
and the quantization procedure for the three-dimensional model of [19] developed in [20]. Here 
we formulate them in a coherent way, which can be applied in quite generic situations. In 
particular, in the next subsection these quantization rules are applied to our model describing 
the degenerate sector of Plebanski theory, and shown to reproduce the correct quantization given 
by the partition function (13.11) . 

We assume that the theory to be quantized has the structure of Plebanski formulation of 
general relativity, i.e. it is represented as topological BF theory supplemented by primary con- 
straints <p, whose time evolution generates secondary constraints ip. Furthermore, we impose the 
partial gauge fixing of the boost gauge freedom, as we did in section 12.21 Then the quantization 
procedure we propose involves the following steps: 

1. First, we need to discretize the primary and secondary second class constraints. We assume 
that their discrete versions give certain restrictions 





(3.20) 
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on the bivectors and the holonomies, respectively. Both constraints are expected to depend 
on the normals x G X assigned to the elements of the two-complex A* and we allowed 
the secondary constraints to depend on the bivectors. In the model we consider here 
this dependence will be absent, which significantly simplifies its spin foam representation. 
However, it is expected to arise in the constraints describing the gravitational sector of 
Plebanski theory. 

2. Using the discrete constraints (13. 20 p . we construct the measures 

V&[B) = A(B,x)6(<f> discr (B,x))dB, 
V^[g] =5^ discr (g,x;B))dg, 

where the first one includes the factor which represents a discretization of the determinant 
of the Dirac matrix, | det {0, Since typically it does not depend on the spin connection, 
it can be expressed through the bivectors and the normals and therefore attributed to the 
first measure only. 

3. Given the measure for holonomies, we evaluate the following quantity, which we interpret 
as vertex amplitude, 

A v (X f ,k ef ,i e ) = I Hv^[g ve ]S {Tv ,x f , ksf ,i e) [g^ {f) g vd{f) ,x e ] . (3.22) 

This is the same formula as (ll.3p given in the introduction, where the boundary state 
is taken to be the projected spin network S(r v ,\ f ,k ef ,i e ) [Sf, x e], and generalizes (13 . 15f) and 
(I3.19p . Due to the gauge invariance of projected spin networks and to the following co- 
variance property of the measure 

V^ B) [g ve g] = P (s ■^ ;sBs " 1) [^ e ], g G Spin(4), (3.23) 

the vertex amplitude (I3.22p is independent of the normals. In contrast, if the measure 
depends on bivectors, this dependence propagates to A v and, as a result, it cannot be 
viewed as a true vertex amplitude in the spin foam representation. If however the measure 
is ^-independent, as it happens in our model, the formula (I3.22p does provide the spin 
foam vertex. 

4. To achieve the correct gluing of the vertex contributions, one has to perform several addi- 
tional steps. The first of them is to evaluate what can be called the vertex amplitude in 
the "connection" representation 

A [g f ,x e ]= [H d ^) II dk °f ) A o( X f» k ef,u)S(p v ,Xf,k ef ,i B ) \gf,x e ] ■ (3.24) 

\f,k ef ,i e XfDv J \{e,f)Dv J 

In generic case this name is not quite precise because this quantity carries a dependence 
on the bivectors, which we did not indicate explicitly, originating from the measure in the 
definition of A v (Xf, k e f, i e ). 
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5. The quantity (13.241) can be already used to glue several vertex contributions. However, 
such gluing can be written more elegantly if one first passes to the ^/-representation. To 
this end, one defines 

A v [B f , x f ]= Y[ [dg/ exp jiTr ((1 + • 9fu(f)gfgja {f )) } MSf, x e ), (3.25) 

J fDv 

where 7 is the Immirzi parameter and we used the fact that the amplitude (I3.25j) depends 
on bivectors only in the combination Bf = 9jdff)^f9fu(j)- The group elements gj e are 
restricted to satisfy the same type of the second class constraints ?/>di SC as g ve . We do not 
integrate over them since they drop out from the final partition function. Alternatively, one 
could insert integrals J T>( x ' B ^{gf e ] which explicitly put the constraints on these holonomies. 
Note in contrast that the group elements gj are integrated with the standard Haar measure 
and not the one involving the second class constraints. 

6. Finally, the total partition function is formed by multiplying the vertex amplitudes (13. 25j) 
using the non-commutative star-productQ [371 EH] and integrating the result over the bivec- 
tors with the measure (13. 21 p 

Z = J U V(x) l B f] (* MBf,x f fj ■ (3.26) 

The order in the non-commutative product is dictated by the orientation of faces and 
by the choice in each face of a "reference tetrahedron" [34J. As above, the covariance 
of the measure ensures the independence of the partition function on the normals, which 
expresses the independence of the quantization on the gauge fixing. 

Several comments concerning this construction are in order: 

• These rules to construct the partition function have been derived by discretizing the canon- 
ical path integral of the original theory and by splitting it into contributions associated 
with different vertices [18] . Therefore, the consistency with the canonical quantization is 
in a sense built-in to this approach. In particular, we will see below how the quantum 
amplitudes introduced above allow to recover various elements of loop quantum gravity. 

• The key element of this construction is the formula for the vertex amplitude (I3.22p . This is 
a straightforward generalization of the usual prescription for the evaluation of the vertex in 
the spin foam models considered in the literature. The only difference is that the measure 
in ( I3.22p is supposed to be non-trivial and, in particular, to include the secondary second 
class constraints. However, this difference has drastic consequences. Without including 
the constraints, the integral defining the vertex is equivalent to the evaluation of a simplex 
boundary state on a flat connection. This recipe has its origin in the topological BF theory 
which describes the dynamics of flat connections. On the other hand, in the presence of 
the secondary constraints in the measure, such correspondence does not work anymore and 
the dynamics becomes more complicated. 



4 This product is denned on the plane waves e 9 (x) = e lTr( - 9 ' x \ with x £ g and momentum given by group 



element g, as e gi * e 
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• While A v (Xf, k e f,i e ) encodes the dynamics, the vertex amplitude in the "connection" rep- 
resentation A v [gf,x e ] describes the kinematical Hilbert space of the theory. It shows that 
the kinematical states always appear as linear combinations of projected spin networks 
[23} IT8] . which is indeed the case for all spin foam models of four- dimensional gravity con- 
sidered in the literature. But not all linear combinations are physically relevant. The space 
of all projected spin networks is too huge and only the states given by (13.241) contribute 
to the path integral. 

• The non-commutative star-product is introduced in the final formula (I3.26P in order to 
combine different exponential factors into the exponential of the discrete BF action. This 
construction is similar to the one introduced recently in the study of a non-commutative 
flux representation [37J |38j [33] . 

• It should be emphasized that at none of the steps one integrates over the normal vectors. 
This is consistent with their interpretation as gauge fixing parameters. Moreover, they 
completely drop out of the partition function, so that in the Lorentzian case an integral 
over these normals would produce an overall infinite factor. 

• An important feature of the proposed quantization is the difference in the roles of the 
primary and secondary constraints. Whereas the latter enter the definition of the vertex 
amplitude and affect the dynamics, the former become relevant only when one glues dif- 
ferent simplices together, i.e. at the very last stage. This is in a drastic contrast with 
the usual spin foam approach where the primary constraints play the central role and the 
secondary constraints have been ignored at all so far. 

3.4 It works! 

Let us now apply the quantization procedure of the previous subsection to our model which, 
after partial gauge fixing, is given by the action (I2.37p . To make it more understandable, we will 
follow the procedure step by step. 

1. The second class constraints in our case are given by and ^ (I2.35p . They have been 
already discretized in section I3TT1 (see (I3.13P ). so that in the chiral notations we have 

<fi disci (B,x) = BJ - XfBjxj 1 , Vdiscr^,^) = (g~ e y 1 x v g^ 1 . (3.27) 

It is useful to note that, besides the second class constraints on the canonical variables, 
the gauge fixed action leads to conditions on the Lagrange multipliers. These conditions, 
provided by (I2.14p and (I2.23p . are naturally combined with the second class constraints, 
which makes possible to extend the latter to spacetime covariant conditions. 

2. The measures (I3.2ip take the form 

vW[B f ] =5(BJ-x f B}x 7 l )dB f , (32g) 

V {x) [g ve ] = 5 {(g^y^vQveK 1 ) d 9ve, 

where we took into account that the factor A(B,x) is trivial due to the last commutation 
relation in ( 12.26p . Note also that the measure over holonomies remains independent of the 
bivectors. 
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3. The vertex amplitude given by the general formula (I3.22p is evaluated in exactly the same 
way as in section [3721 Similarly to (I3.19p . one finds 



A v (X f ,k ef ,i e ) = {^]}X\d k H ku{f)fkd(f)r 



(3.29) 



Thus, up to a normalization factor, we again reproduce the correct vertex amplitude of 
the SU(2) Crane- Yetter model fl3~TD . 

4. The rest of the construction will restore the correct face and edge amplitudes. But not 
only that. It will also teach us several important lessons. In particular, let us evaluate the 
vertex amplitude in the "connection" representation C I3. 24D . Since the vertex (13. 29ft does 
not depend on A/, the sum over this label can be done explicitly. Indeed, it is easy to 
prove the following identity 



m,m' ,n,n' 

due to which one finds 



n4 / 5(g / -(x)( g ;(a;))- 1 ) 



(3.31) 



where we introduced g/(x) = gx^ (f) Sf9x d(f) and S(r v ,k f ,i e ) is the usual SU(2) spin network 
associated with the boundary graph IV This result has two remarkable features. First, 
the delta function on the r.h.s. imposes the same condition as the second class constraint 
(I3.27P on the holonomy g u (f)d(f) • Second, the amplitude (13.311) is represented as a sum over 
SU(2) spin networks, which shows that the latter span the kinematical Hilbert space of 
our theory in full agreement with the SU(2) Crane- Yetter model. From this one concludes 
that the sum over the auxiliary representation labels, which the vertex amplitude does 
not depend on, has two important effects: on the one hand, it restores the secondary 
second class constraints on the arguments of the boundary states and, on the other hand, 
it provides the reduction from the space of all projected spin networks to the kinematical 
Hilbert space of the theory. We emphasize that this reduction does not involve the primary 
simplicity constraints at all and is achieved due to the secondary constraints encoded in 
the form of the vertex amplitude! 



5. The next step is to substitute (13.31 j) into (I3.25p . This gives 



A v [B / ,x / ] = Y f {mU 

kf,i e fDv 



SU(2) 



S(r v ,kf,i e )[hf), 



(3.32) 



where we used the chiral decomposition of the 50 (4) trace ( 1A.7|) . the gauge invariance of 
spin networks, and denoted 



«u(2) 9 bf = \ (1 + t- 1 )^ + \ 



1 - ; )x,'//,x ; . 



(3.33) 
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6. Finally, we glue the amplitudes (I3.32p together by means of the formula (I3.26p . To distin- 
guish the representations and the group elements associated to different vertices, we put 
the index v on them, i.e. in (I3.32p one should make the replacements (kf,i e ) >->■ (k v f,i ve ) 
and h j i — y h v f. Then the non-commutative star-product ensures that each face comes with 

the factor exp |itr (b^Hf)\ where Hf = Ylvcf^vf- Since the geometric meaning of h v f 

is the positive chiral part of the curvature around the wedge (vf) (see Fig. CD), the group 
element Hf gives (the positive chiral part of) the full curvature around the face /. On 
the other hand, due to the primary simplicity constraints entering the measure on the 
bivectors, one has = so that the integrals in (13.261) generate S(Hf) imposing the 
flatness condition. Expanding the 5-function in the sum over representations, it is easy to 
see that the full partition function is given by 



Z 



n 

jf k v f,i ve {v,f) 



d kvf I dh v f 

SU(2) 



n 

/ 



d 3fXj f ( n 

.vCf 



] [ [{ 15 J} s (T v ,k vf ,i e )[h vf 



(3.34) 

where Xj is the SU(2) character of representation j. It is immediate to check that, doing the 
remaining integration and contracting all indices, one reproduces the SU(2) Crane- Yetter 
state sum (13.11) with the same face, edge and vertex amplitudes. 



Thus, we conclude that the quantization rules given above lead to the correct spin foam 
quantization of the degenerate sector of Plebanski theory and therefore provide the correct 
implementation of all the constraints. 



3.5 Why does it work? 

In the previous subsection we went through a long way to get the partition function of the 
constrained theory. In fact, in the particular case of our model, there is a shorter way to arrive 
at the same result, which partially explains the origin of the proposed quantization rules. It relies 
on the observation that the vertex amplitude in the "connection" representation A w [g/,x e ] has a 
much simpler expression. Indeed, in [18] it has been proven that its spin foam like representation 
fl3T2^j) follows from 

MSf,^} = f U V(X '' B) ^e} U 6 0W)g/^(/)) • (3 - 35) 
eDv fDv 

As a result, one immediately obtains 

A v [B f ,x f } = J n^ (z;B) [^e]n ex P( iTr [( 1 + ^)5/-G', / ]), (3.36) 

eDv fDv 

where G v f = g fu(f)gu(f)vgvd(f)gd(f)f is the curvature around the wedge (vf). Then the partition 
function (I3.26P produces the exponential of the unconstrained Spin(4) BF action integrated 
with the measure (13.211) involving the primary and secondary second class constraints. This 
is nothing else but a discretization of the standard canonical path integral for the constrained 
theory. In our case the discrete constraints (I3.27P can be solved explicitly, which simply reduces 
the path integral to the SU(2) sector, and therefore the coincidence with the reduced phase space 
quantization (13.11) is guaranteed. 
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One can ask: why did we do all the above complicated calculations if they are not required 
to get the final spin foam model? The point is that the shorter way is available only if it is 
possible to explicitly find the reduced phase space at the discrete level. Although this can be 
done for our simple model, this seems to be out of reach in more complicated situations such 
as the gravitational sector of Plebanski theory. On the other hand, the representation ( 13.24ft 
disentangles the vertex contributions and the gluing of different vertices, and importantly this 
is done before implementing the constraints. This can be viewed as the first step towards the 
spin foam representation of the partition function, which should follow after integrating out the 
remaining geometric variables. 

Furthermore, the derivation of the previous subsection clarified many subtle issues such as 
the imposition of constraints, gauge invariance, the role of the Immirzi parameter, etc. Some of 
them have been already discussed above, and we will summarize once more our main conclusions 
and observations in the next section. 

4 Discussion 

In this paper we studied the classical and quantum descriptions of the degenerate sector of 
Plebanski formulation of general relativity. We have shown that one of its subsectors, analogous 
to the gravitational sector of Plebanski theory, provides an interesting and useful model to test 
the ideas of the spin foam quantization. Classically, it represents a constrained four- dimensional 
Spin(4) BF theory which, upon elimination of the constraints, reduces to the SU(2) BF theory. 
Since both these theories are of BF type, their spin foam quantization is well known and can 
be used to find the correct way of implementing the constraints. As a result, we formulated a 
general procedure to build the partition function of any constrained theory of Plebanski type. 
In particular, we showed that, being applied to our model, this procedure works perfectly, giving 
rise to the right kinematical Hilbert space and generating the right dynamics, i.e. those which 
agree with the SU(2) Crane- Yetter model. 

One of the main results of this analysis is the clarification of the role of the primary and 
secondary second class constraints in the construction of the spin foam partition function. As 
has been already argued before [18J [16j [19j [20], the secondary constraints affect the measure for 
holonomy variables and determine the form of the vertex amplitudejf] On the other hand, the 
primary constraints enter only at the very last step of the construction when the contributions 
of different simplices are glued together. In fact, this is a very natural result. The primary 
constraints is a simple consequence of the choice of our basic variables Bf, which live on the 
boundary of 4-simplices. At the same time, the secondary constraints appear as a commutator of 
the primary ones with the Hamiltonian and therefore contain information about the dynamics of 
the theory. Moreover, they constrain the variables g ve living "inside" 4-simplices. Giving these 
observations, it should not be surprising that the main quantity responsible for the dynamics 
in the spin foam approach is governed by the secondary constraints, whereas the primary ones 
play only some minor role at the boundary. Furthermore, even the kinematical Hilbert space 
of the constrained theory turns out to be completely determined by the secondary constraints 
because the kinematical boundary states appear as projected spin networks weighted by the 
vertex amplitudes (see ( I3.24p ). 

All these statements are in a drastic contrast with the usual constructions performed in 

5 A similar modification of the measure for holonomies has been found also in recent group field theory 
constructions [34] . 
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the four-dimensional spin foam models of general relativity where the principal role is given to 
the primary constraints, whereas the secondary constraints are not considered at all. In the 
quasiclassical limit these models do provide the right dynamics since in this limit one sets on 
shell where the secondary constraints are effectively induced. However, beyond the limit the 
primary constraints are not sufficient to suppress the quantum fluctuations of the degrees of 
freedom constrained by the secondary ones. This is why the latter should be taken into account 
and are crucial to get the right dynamics at quantum level. This is clearly demonstrated by our 
model as well as by its three-dimensional analogue [19]. 

It is worth also to note that our derivation confirmed once more that the closure constraint 
of Regge calculus should not be imposed in spin foam models. This constraint requires the 
invariance of intertwiners which is achieved by integration over the normal vectors x e associated 
to tetrahedra of the simplicial decomposition. However, our model, in agreement with previous 
claims [23] IT8] [39] . clearly shows that such integration would be inconsistent and the normal 
vectors should be kept fixed, which is nothing else but the usual gauge fixing of the boost gauge 
freedom in the gravity path integral. 

Another interesting point is that our quantization has been done in the presence of the 
Immirzi parameter. Nevertheless, it did not have any effect on the constrained theory both at 
classical and quantum level. After imposition of the constraints, it completely drops out of the 
partition function. This should be compared with the claims that its appearance in loop quantum 
gravity is a consequence of an unfortunate choice of variables (the Ashtekar-Barbero connection 
which is not a pull-back of a spacetime connection) and with a right choice it disappears from 
physical results [UJ [2SJ [T7]. 

Finally, let us comment on the extension of our construction to the physical case of the 
gravitational sector of Plebanski theory. The main difference distinguishing it from our model is 
the form of the constraints which cannot be written anymore in the simple form f 1 2 . 3 5 1) or ( I3.27P 
after discretization. In particular, the secondary constraints become explicitly dependent on the 
.B-field [HI H2] . Although the construction of section 13.31 is still well defined in the presence 
of such dependence, it gives rise to many complications. The most important one is that the 
quantity ( 13 .22 j) starts to depend on the bivectors Bf and its interpretation as a vertex amplitude 
is not viable anymore. It is not clear whether this is a serious problem or just a minor obstacle. 
In principle, the 5-dependence can make impossible to integrate out the bivectors because the 
resulting integrals are not of BF type anymore. However, given that in our procedure this 
integration appears only as a way to glue the simplex contributions, one may still hope that it 
will be possible to factorize the -B-dependence and to extract the spin foam vertex. In any case, 
this issue certainly deserves a further study. 
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A Conventions 



A.l so (4) algebra and the chiral decomposition 

Our conventions for indices are the following: /i, v, . . . denote spacetime indices, a,b, . . . are 
spatial indices, I, J, . . . are the so(4) indices in the tangent space, and i, j, . . . refer to the su(2) 
subalgebra. e 1JKL aric [ e l i k are the invariant fully anti-symmetric tensors in four and three 
dimensions, respectively, normalized by e 0123 = e 123 = 1. Besides, we use (••) and [••] to denote 
symmetrization and anti-symmetrization, respectively, with weight 1/2. Since all tangent space 
indices are contracted with a unit metric, we do not follow the rule, that the contracted indices 
should be in opposite positions, very strictly. 

The canonical basis of so (4) is composed of the rotation and boost generators, Lj and Lfj, 
which have the following commutation relations 

[Li, Lj] = Eij L k , [Ki, Kj] = Eij k L k , [Ki, Lj] = Eij k K k . (A.l) 

The rotation generators form the canonically embedded su(2) subalgebra. On the other hand, 
given a four- dimensional normal vector x 1 , one can introduce a boosted subalgebra su(2) x . It 
is formed by the generators which leave the vector invariant. If one introduces the covariant 
notation for the so (4) generators T lj = \e^ k L k , T° l = \K l , then 

/£(*) = ~ e ijmp e klnp x m x n = 5 J K J L - 2^4^], 
J I K J L (x)=2x^ J S I [ ] K x L] . 

are the projectors on su(2) x and its orthogonal completion, respectively, i.e. acting on the 
generators in the vector representation, they satisfy Ij/ L (x)T KL ■ x N = 0. 

The antisymmetric bivectors B IJ form the adjoint representation of the so (4) algebra. On 
this representation we define the action of the Hodge operator as {*B) IJ = \ e IJ klB kl . Since 
the Hodge operator squares to one, * 2 = id, it splits the space of bivectors into the direct sum 
of two eigenspaces with eigenvalues ±1, 

B = B\ +) Tl +) + B { r ] T t { -\ *rW = ±rW. (A.3) 

This corresponds to the chiral decomposition Spin(4) = SU(2) x SU(2), and the generators of 
the Lie algebras of the two chiral subgroups are given by = \ (Lj ± Ki) with 

Pf \ ?f } ] = % fc 7f } , [Ti + \Tj~^] = 0. (A.4) 
In the adjoint representation we normalize the trace such that 

Tt(T ij T kl ) = 5 ! K J L =j> Tr(AB) = A IJ B U , (A.5) 
where = S^-S^j. Then the chiral decomposition ( 1A.3I) implies that 

Tr(^ (e) Tf } ) = 6^, =^ Tr(AB) = A { + ] B\ +) + A^ B\~\ (A.6) 

On the other hand, it is useful to remember that the su(2) generators Lj have a different nor- 
malization, namely, tr (LiLj) = 25y. Due to this, in terms of the su(2) elements B^' = B^L 1 , 
the so (4) trace reads 

Tt(AB) = X - tr (A i+) B {+) ) + - tr (A H 5 H ). (A.7) 
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A. 2 Clebsch-Gordan coefficients 

Our conventions for SU(2) invariant intertwiners follow [H]. A generic invariant intertwiner is 
denoted by i and is supposed to be normalized as 

mi-mi 

In the particular case of three coupled representations, the matrix elements of the intertwiners 
are given by the Clebsch-Gordan coefficients C 3 ! 33 . It is convenient also to define the two 
invariant maps based on these intertwiners 



C nnn : H Ul) <g> H U2) -> H ij3) and C nnn : U^ z) ^H Ul) ®U {j2) '. (A.9) 

Then the following properties are satisfied 

h+h 

C ^C 31323 ' = dj 1 5 jf l„ Y. d 3 C 31323 C 31323 = l n ®l n , (A.10) 

3=\h-h\ 

where dj — 2j + 1 is the dimension of the SU(2) representation. In particular, the first property 
ensures that one gets the right normalization (1A.8j) . Finally, we fix the normalization of Wigner 
matrices by requiring 



/ dhD$ n (h)D%l(h) = dfSjfS^S^, (A.n; 

JSV(2) 



'SU(2) 

With these normalizations, the matrix elements are recoupled as follows 

oa»DS 2 (/ l )= x! Y, d AtL<t>$ n {h). (A.i2) 

3=\h-h\ m ' n 

These properties are sufficient to prove (13.181) and (13.301) . 



B Canonical analysis of the degenerate Plebanski sector 

The canonical analysis of the gravitational sector of Spin(4) Plebanski theory has been carried 
out for the first time in and further elaborated in [321 S3]- The degenerate sector described 
by the action (12.11) can be analyzed along a similar way. Here we follow the original method of 

EH- 

One starts as usual from the 3+1 decomposition. One immediately observes that u^ J and 
Pjj = e abc Bfc appear as conjugate variables, whereas the action does not contain time derivatives 
of A MiypfT , 0Jq 3 and Bq 3 which therefore are expected to play the role of Lagrange multipliers. How- 
ever, the latter variables appear quadratically in the term generating the simplicity constraints 
which leads to certain complications in considering them directly as Lagrange multipliers. To 
avoid these complications, one adds new non-dynamical variables /i^" 7 and njj which enforce the 
vanishing of the momenta conjugate to by means of the following additional term 

J^d\(Tr(n a d B 0a ) - Tr(/i a 7r a )) . (B.l) 
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Then the pase space is spanned by u^ J , and 7r^ J with the symplectic structure given by 

K J (x), P b KL (j)} = ^0(x,y), {B^y^iy)} = S^J^y), (B.2) 
and the total action leads to the following primary constraints 



$(B,B) ab = -£ IJKL B^BQ b 

msr, = \e IJ KL P?X b 





(B.3a) 




(B.3b) 


o, 


(B.3c) 


0. 


(B.3d) 


0. 


(B.3e) 



The last three are nothing else but the various components of the simplicity constraints (12. 2p . 
But as is shown in appendix [Cj the six constraints §(B,B) ab are in fact linear combinations 
of $(P, P) a b and $(P, P) ab , and do not require a special attention. The Hamiltonian can be 
written in the following form 

-H = e ahc B^Fl J + ul J g u + A 0a06 $(P, B) ab + \ X 0abc e bc MP, B) d a 

2 (B.4) 

+ Yq \ abcd e abf e cdg $(P, P) f9 ~ »a*h- 

The next step is to study the conditions imposed by the conservation of the primary con- 
straints. Let us compute their time derivatives by commuting them with the Hamiltonian (1B.4[) . 
First, one finds 

TTjj = s ahc Fl 3 c + 2A 0a0b (*P)^ J + X - \ 0abc e bcd ^P)jj « 0. (B.5) 

These 18 equations split into two sets. By contracting them with P}j, the last terms produce 
the simplicity constraints so that one remains with 9 conditions 

C ab = ^cdpljpbj w Qj (B.6) 

which should be interpreted as secondary constraints. The remaining 9 equations following from 
flB.51) can be written as 



e bcd Ti{P a * F cd ) + 2A 0Wc Tr(P a P 0c ) + - X obcd 6 cdg Ti{P a P g ) = (B.7) 



and fix the Lagrange multipliers A 06cd . 

The conservation of the Gauss constraint Qu does not generate new conditions provided it 
is shifted as follows 

Gij = Gij+[B 0a ,n a ]jj. (B.8) 

The shift ensures that Q' u is a generator of Spin(4) gauge transformations and since the Hamil- 
tonian is gauge invariant, Q'jj is preserved under evolution. 

Next it is convenient to consider $(P, P) ab . Its commutator with the Hamiltonian generates 
6 new conditions 

^ab = £ cd(a Ti ^ Bod ^ Dc pb)j w o, (B.9) 
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which give rise to additional secondary constraints. Having obtained these constraints, we can 
now prove a very useful Lemma which facilitates a lot the following analysis. 

Lemma: Let b 1 be defined by the solution of the simplicity constraints in any of the two 
sectors (see ( 12. 3 p ). Assume that i) p^ = e abc bl c is invertible in the sense that there exists p a such 
that 

PjP b i = Sl p a J = 5 I J -x I xj- (B.10) 

ii) the matrix 

Q ab ' cd = e IJK L x I p ( je b)9{c p l ^ g (B.ll) 

in invertible. Then D a x can be expressed as a linear combination of \P ' , Q a = Tr(P a -kQ) and 
the simplicity constraints, and therefore it weakly vanishes. 

Proof: First, it is easy to see that in both solution sectors (I2.3p . i.e. on the surface of the 
simplicity constraints, one has 

Q a =~e IJKL x I p a J p b K D b x L , 

2 (B.12) 

qab = e^Eu^X^iffilD^L. 
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Then one can check that 



Q^plD^X! = 2V ab - g {a s b ^ d bi c p dJ . (B.13) 

Since the matrix Q ab > cd is assumed to be invertible, this implies that p^DyXi is weakly vanishing. 
Finally, one verifies that 

D a xi = p) {p J {a D b) xj) - Eijklx j pJpjQ h « 0. (B.14) 

□ 

We still need to analyze the stability of two simplicity constraints, ( IB.3ch and (IB.3dh . How- 
ever, since the former is expressible through the latter and (lB.3ep . only $(P, B)\ remains to be 
considered. Its time derivative leads to the following condition 

$(P,B) a b = -Tr([uj ,B ob }*P a ) + Tr(ii b *P a ) « 0, (B.15) 

where we neg lected the term 2e acd Tr(B 0b * D c B od ) because, on the surface of the simplicity 
constraints, it is proportional to D c x T and vanishes weakly by the above Lemma. The resulting 
equation fixes 9 of the 18 components of the Lagrange multiplier /i^ J and does not lead to new 
constraints. 

This completes the analysis of the primary constraints. But the appearance of the secondary 
constraints C ab and ^ ab requires to repeat the stabilization procedure. However, at this step it 
works rather differently for the deg-gravitational and deg-topological sectors. Due to this reason, 
we consider them separately. 

B.l Degenerate gravitational sector 

First, let us consider the conservation of the secondary constraints ^ ab . This gives the following 
equation 

qab = 2£ «*(» £ ^Tr(£W * D c D f B 0g ) - £ cd(a Tr([ Wo , B od ] * D c pV) + ^Tr(^ * D c pV) 
+ \ 0cd ( a Tr(pV[B Qc , B Qd \) + ~ \^E cfg E rpq E cd{a i:T{P b \P\B M }) « 0. 

o 
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However, it is easy to see that the first term weakly vanishes due to the Lemma, whereas the next 
two terms produce the equation flB.15j) plus contributions proportional to D c x l . As a result, the 
stability condition reduces to the vanishing of the last two terms. The crucial question for us is 
the form of the matrix in front of the Lagrange multiplier A^ 9P<? in the last term, which arises 
from the commutator of ^/ ab with the primary constraints $(P, P) cd . In the deg-gravitational 
sector, where the P-field is given by ( 12.3aj) . it is found to be 

{^! a \ $(P, P) cd } = AP'fj e b)9{c P} ] K B^ g S v Q ab ' cd . (B.17) 

Since this is the same matrix which appears in the Lemma and generically it is invertible, 
the condition flB.16j) fixes the six Lagrange multipliers X f9pq . At the same time, it shows that 
$(P, P) ab and ^ ab are mutually non-commuting. 

Before we proceed further, we prove an additional useful result that the curvature of the spin 
connection is weakly vanishing. Indeed, using fl2.3al) and the inverse field introduced above, 
one has 

£ a b c F IJ = ^cjIJ^pKL + ^cjU^KL 

= e lJ KL x K y L h C ab + 4x [J e abc D b D c x I] « 0, 

where in the second term we represented the curvature as a commutator of two covariant deriva- 
tives. Evaluating now the time derivative of the secondary constraint C ab 

C ab = 2s acd £ b ^Tr(B 0g D f F cd ) - 2A°^ c Tr(P b * D c B og ) - l - X^e cfg e rpq e cda Ti(P b * D d P r ) « 0, 

(B.19) 

one immediately concludes that the stability condition is satisfied due to the above Lemma and 
the vanishing of the curvature proven in ( IB. 181) . 

As a result, the stabilization procedure finishes at this point and the list of all constraints 
is given by nfj, Q u , ®(P,B) a b , $(P, P) ab , ^ ab , and C ab . Note however that the last constraints 
are reducible. Namely, due to the Bianchi identity, they satisfy 

D a (C ab pi) = I e IJKL e acd F c K d L D a xj, (B.20) 

where the r.h.s., as we know, vanishes on the surface of the other constraints. Thus, only six 
components of C ab are independent. To split the resulting constraints into first and second class, 
we introduce 

tt$ = Tr(vr a P b ), 7rg = Tr(7r a * P b ). (B.21) 

Then, using the weak vanishing of the curvature and of the covariant derivative of the normal 
vector, it is straightforward to verify that 7r^, Q\j and C ab are first class, whereas tt^, $(P, B) a b , 
$(P, P) ab and ^ ab are second class. This implies the following counting of degrees of freedom. 
The original phase space is 4 x 18 = 72 dimensional. The second class constraints remove 
9 + 9 + 6 + 6 = 30 degrees of freedom, whereas the first class constraints together with the 
corresponding gauge fixing conditions fix 2 x (9 + 6 + (9 — 3)) = 42 of them. This leaves us with 
a zero-dimensional phase space confirming that we are describing a topological theory. 

If one partially fixes the gauge, taking the normal vector x 1 to be a given function x 7 (x), this 
gauge fixing condition can be combined with $(P, P) ab to get the simplicity constraints (I2.13bl) . 
Moreover, the part of the Gauss constraint Q a is conjugate to the gauge fixing condition and 
therefore becomes second class. As shown by our Lemma, it can be combined with ^> ab to get 
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the secondary constraints (12.241) . Besides, the constraints &(P,B) a b become identical to the 
condition f |2.14j) on the Lagrange multipliers, whereas C ab and the remaining part of Qu are 
equivalent to Cf and Qi defined in (12.1 7[) and (12.281) . respectively. After integrating out the 
auxiliary variables iijj, one recovers the same phase space and the same constraint structure as 
in section l2~2l where the canonical analysis was performed for the gauge fixed action (I2.10p . 
Furthermore, it is possible to show the following equality 

| det {x, G}\ 5(n (1) )5(n {2) )5 ($(P, P)) 5 ($(P, B)) 8{*)8{Q)8{x - x(x)) 
~ 8Wj)8(xjB£)8(x?P!j)6{DJ), 

where 

A D ~ (detp) 6 (detg) 2 (B.23) 

is the determinant of the Dirac matrix of the commutators of the second class constraints, the 
second factor is a part of the Faddeev-Popov determinant corresponding to the gauge fixing of 
the normal x 1 , and all equations are given up to numerical factors. The meaning of this result is 
twofold: first, it demonstrates explicitly the recombination of the constraints mentioned above 
and, second, it confirms the triviality of the factor A from (I3.2ip (see (|3.28|) ). Thus, all results 
derived from the gauge fixed action (I2.10p are in the full agreement with the complete canonical 
analysis of the initial gauge invariant theory. 



B.2 Degenerate topological sector 



Now we turn to the deg-topological sector described by the solution ( I2.3b|) . The conservation 



of the secondary constraints ^ ab still generates the condition (IB.16j) where the first line can 



be dropped. However, in contrast to the previous case, the matrix appearing in the last term 
vanishes 

{ty a \ $(P, P) cd } = pf" e b)9{c Pf K Blf ^ P (B.24) 

so that $(P, P) an d ^ are now mutually commuting. Moreover, the forth term also vanishes 
in the deg-topological sector so that the stability of ^ ab does not generate any new conditions. 
The second crucial difference arising in this sector is that the secondary constraints C ab ( 1B.6[) 



turn out to be linearly dependent with other constraints. Indeed, using ( I2.3bj) and the same 
trick as in (1B.18|) . one obtains 

C ab = -2e? cd $D c D d x I . (B.25) 

Then the above Lemma ensures that this quantity can be expressed as a linear combination of the 
Gauss constraint, simplicity and \E ,ab . Due to this, C ab does not require a separate consideration 
and its stability follows from the above analysis. 

As a result, the independent set of constraints is provided by Tijj, Qjj, $(P, B) a b , $(P, P) ab , 
and ty ab . To split them into first and second class, besides (IB . 2 1 1) . we define 

yab = ^afe _ n cd} (^(PP)" 1 ) $(P, B) f c . (B.26) 

Then one can check that tt^), Q'iji $(P, P) a& and ^ ab are first class, whereas tt^ and $(P, B) a b 
are mutually non-commuting and therefore give second class constraints. The counting of degrees 
of freedom works as above. The second class constraints give 9 + 9 = 18 conditions, whereas 
the first class constraints and their gauge fixing conditions produce 2 x (9 + 6 + 6 + 6) =54 
more. Altogether they fix all degrees of freedom of the original 72-dimensional phase space, also 
showing that this sector describes a topological theory. 
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C Constraints for constraints 

Let us consider the simplicity constraints split according to the 3+1 decomposition as in (IB.3|) 
and written in terms of the chiral variables. They take the following form 

HB,B) ab = B^B^ - Bt/B^\ 

HP, B)\ = Pl +)a B^ - Pt )a B^\ (C.l) 
$(P P) ab = p(+) a p(+) b _ p(-) a p(-) b 

Then it is straightforward to check that the following linear combination of constraints 

\ (~ Pfr^Pt^iB, B) ab - B£*P£**{P, Bfl 

-p£ )ci B^ h m Bfi + b£*b<-**{P, p>*) p^p^p^p^ ( c -2) 

p(+)»a p(+)c2 p(+)c 3 R (-)u p(-)da pH**/p p\ci d i 
"-°0a ^i! r i 3 -°06 ^2 r is ^\ r i r l 



ab & c\C2C-^ d\d,2<lz 



identically vanishes. This can be done either expanding explicitly all sums over repeated indices 
or assuming that P\ are two invertible matrices which allows to write, for instance, 

e clC2C3 P/ 2 +)c2 P/ 3 +)c3 = ( det P«) £il -3 (P( + H); . (C.3) 

Using this property in each term in (1C.2I) . the vanishing of T a f, follows trivially. 

The nine quantities T ab encode constraints for the simplicity constraints f)2.2p . However, 
they themselves are not all independent. If one considers their antisymmetric part Tr^i, it is 
possible to show that the coefficients in front of the constraints vanish on the constraint surface. 
Indeed, assuming again the invertibility of Pj , one finds 

6T H = (PW- 1 ); - (P^-%) ($(P, B f a] - B ™ (P^-% m PY d ) , (C4) 

where the prefactor can be rewritten as a linear combination of the original constraints ( 1C.1I) 

2(tr(p(+)p<+))- 1 )^(<l>(P,B) i 6 -B<-) i (pW- 1 )-*(P,P)»"'). 

Thus, only the symmetric part of ( 1C.2I) generates constraints for constraints and allows to express 
§(B,B) a b in terms of other 15 simplicity constraints. 

One can also notice that the vanishing of ( 1C.2I) continues to hold even if one replaces there 
one or two of the fields B 0a by P a . For example, replacing only one field, one obtains the 
following combination of constraints 



T a — r p 

1 b c ciC2C3^did,2d3 



1 Q M Cldl $(P, B)\ - M adl $(P, B) c \ - P^B^^P, P) adl 



2 



+P l i+)a B { ob )l <S>(P,P) cldl ) M C2d2 M c:id3 - P l i+)c2 B { ~ )i M ad2 M c:ida ^(P,P) Cldl 



(C.6) 



where we denoted M ab = ~ tr (p(+)p(~)) a6 . This combination is also vanishing, so that one 
could think that these are additional constraints for constraints. However, in contrast to (IC.2I) . 
T a b vanish for any $(P, B) a b and symmetric $(P, P) ab , not necessarily of the form (IC.lj) . Again 
this can be easily proven by assuming the invertibility of the matrix M ab . Due to this, they do 
not reduce the number of independent constraints which remains to be 15. 
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